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We discuss the flow past a flat heterogeneous solid surface decorated by slipping stripes. 
The spatially varying slip length, b(y), is assumed to be small compared to the scale of 
the heterogeneities, but finite. For such weakly slipping surfaces, earlier analyses have 
predicted that the effective slip length is simply given by the surface-averaged slip length, 
which implies that the effective slip-length tensor becomes isotropic. Here we show that 
a different scenario is expected if the local slip length has step-like jumps at the edges 
of slipping heterogeneities. In this case, the next-to-leading term in an expansion of the 
effective slip-length tensor in powers of max(6(y)) becomes comparable to the leading- 
order term, but anisotropic, even at very small b(y). This leads to an anisotropy of the 
effective slip even for very weakly slipping stripes, and to a significant reduction compared 
to the surface-averaged value. The asymptotic formulae are tested by numerical solutions 
and dissipative particle dynamics simulations. 



1. Introduction 

With emerging technologies in microfluidics ( Stone et aZ.ll2004 ; Squires fc Qua ke 2005), 
there has been renewed interest in quantifying the effects of surface chemical hetero- 
geneities with local scalar slip ( Vinogradova! 19991 Lauga et al. 2007 ) on fluid motion. 
Well-known examples of such heterogeneous systems include composite su perhydropho - 
bic (Cassie) surfaces, where a gas layer is stabilized by a rough wall texture (|Querell2008f h 
For these surfaces effective slip lengths are often very large comp ared with the value on 
flat solids (jBocquet fc Barrat" 20071 ; Vinogradova fc Belvaev 2011 ). so that they have the 
potential to influence microfluidics (or to extend microfluidic sy stems to nanofluidics), 
by ge nerating very fast and well-controlled flows in small devices ( Vinogradova fc Dubovl 
20121) . 

In such situations it is convenient to construct so-called effective slip boundary con- 
ditions, where the complex flow pattern at a heterogeneous surface is repl aced by an effec- 

five flow averaged oyer th e length scale of the experimental configuration ([Vinogradova fc Belvaev 
2011 ; Kamrin et a/.ll2010h . In other words, rather than trying to solve equations of motion 
on the scale of the individual corrugation or pattern, one considers the macroscale fluid 
motion (on the scale larger than the pattern characteristic length) by using macroscopi- 
cally equivalent boundary conditions for an imaginary smooth surface. Such an effective 
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Figure 1. Sketch of the striped surface: = n/2 corresponds to transverse stripes, = to 

longitudinal stripes. 



condition mimics the actual one along the true heterogeneous surface. It fully character- 
izes the flow at the real surface and can be used to solve complex hydrodynamic problems 
with much reduced computational effort. The effective slip appr oach has been supported 
by statistical diffusion arguments ( Bazant fe Vinogradova 20081) . justified for the case of 



20101). and recently confi rmed 



Stokes flow over a broad class of surfaces jKamrin et all 1 _ . 

by computer simulations (jPriezievll2011 ISchmieschek et aZ.ll20ia IZhou et aLll2012h . 

For an anisotropic texture, the effective boundary condition generally depends on the 
direction of the flow and is a tensor, b e fi = {bff}, represented by a symmetric, positive 
definite 2x2 matrix, which can be diagonalizcd by a rotation with angle (Figure [T]). For 

all anisotropic surfaces its eigenvalues &jj ff and b^ s correspond to the fastest (greatest for- 
ward slip) and slowest (least forward slip) orthogonal directions (IBazant fe Vinogradova 
2008). In the general case of arbitrary direction 9, the flow past such surfaces with 
anisotropic effective slip becomes misaligned with the driving force. Therefor e, anisotropic 
textures can potentially be used to generate transv erse hydrodynamic flow (IAidarill2002 ; 
Bazant fe Vinogradovall2~008t Stroock et aZ.|[2002a| ) , which is of obvious fundamental and 
practical interest. For example, transverse hydrodynamic couplings in flo w through a tex- 
tured channel can be used to separ ate/ concentrate suspe nded particles (|Gao fe Gilchrist 
2008 ) or for passive chaotic mixing ( Stroock et "aZll2002alf6h . This can also be used to gen- 



erate anisotr opic electrokinetic flows (jBahga et al\ l2010t iBelvaev fe Vinogradova! 12011 
Sauiresll2008h . 



However, it has been predicted that regardless of the anisotropy of the surface texture, 
the effective slip-length tensor, b e fi, becomes isotropic (b^ s = &jj ff ) for a weakly slipping 
pattern, i.e. when the local slip length, 6, is small compared to the characteristic scale of 
heterogeneities, L. The value of the effective slip is the surface average of the local slip 
length, b c ff = I{b(x,y)}. In the particular case of a no-slip plane covered by patterns 
with constant slip length b — the s ituation considere d in most previous pub lications on 
the subject (IPrieziev et all 2005 ; Ybert et al. 2007 ; Feuillebois et al. 20091 ) - one can 
easily derive (jBelvaev fc Vinogradova 2010a ; Kamrin et all 2010 ) 
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(i.i) 



where — S/L is the surface fraction of the slipping phase. This result implies, among 
other, that the flow aligns with the applied driving force fo r all in-plane directions. Thus 
it seems impossib le to generate transverse hydro dynamic ([Vinogradova fc Belvaev 2011 



Zhou et all 



20121 ) or transverse electro-osmotic ( Belvaev fc Vinogradova! 2011 ) phenom- 



ena for weakly slipping anisotropic textures. Another important, and somewhat remark- 
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able, consequence of Eq. is that the effective slip is predicted to depend only on 

the fractions of slipping areas, but not on their detailed structure. Known derivations of 
Eq. however, neglect localized flow perturbations around possible jumps in discrete 

slip lengths, from to b, at the border of heterogeneities. 

In the present paper we reconsider this problem, focussing on the situation where 
the perturbation b(x,y) is piecewise constant, i.e., it jumps in a step-like fashion at 
heterogeneity boundaries. In Sec. [5] we define the problem and pro ve that known two- 
term expansions for effective slip lengths of one-dimensional textures ( Kamrin et a/.ll2010l) 
are not valid for discontinuous local slip. Then in Sec. [3] we construct the expansions for 
the eigenvalues of the slip-length tensor of alternating weakly slipping stripes. Here we 
also analyze a singularity of the velocity gradient at the edges of stripes. The details of 
the computer simulation method (dissipative particle dynamics) are discussed in Sec. HI 
Finally, in Sec.0 we present simulation and numerical results to validate the predictions 
of the asymptotic theory. The practical implications and limitations of our models are 
also reviewed here. A brief description of the numerical scheme is given in Appendix [A] 



2. Problem set-up 

We consider a creeping flow along a planar anisotropic wall, and a Cartesian coordinate 
system (x, y, z) (Figure [TJ. The origin of coordinates is placed at the flat interface, a ID 
texture varies over a period L. Our analysis is based on the limit of a thick channel or 
a single interface, so that the velocity profile sufficiently far above the surface may be 
considered as a linear shear flow. Dimensionless variables are defined in terms of the 
reference length scale L, the asymptotic shear rate far above the surface, G, and the fluid 
kinematic viscosity, v. 

For a one-dimensional texture there exists a simple relation between longitudinal and 
transverse effective slip lengths ( Asmolov k, Vinogradoval 2012 ) 



_4l£M. ,2,) 

Therefore, it is sufficient to consider the longitudinal configuration. Since in this case, 
the velocity has only one component, we seek a solution for the velocity profile of the 
form 

v = U + u, 

where U = z is the undisturbed linear shear flow. The perturbation of the flow u (y, z), 
which is caused by the presence of the texture and decays far from the surface at small 
Reynolds number Re = GL 2 /v, satisfies the dimensionless Laplace equation, 

Au = 0. (2.2) 

The boundary conditions at the wall and at infinity arc defined as 

z = : u-ef3(y)d z u = e/3 (y) , (2.3) 
z ->• oo : d z u = 0, (2.4) 

where e = b max /L and j3 — b (y) /6 max is the normalized slip length. 

For a weakly slipping anisotropic texture, £ < 1, an expansion in powers of e yield s 



to leading order an area-averaged isotropic slip length (jBelvaev fc Vinogradovall2010al ). 
In practice, this means that the slip-length tensor becomes isotropic and that for all 
in-plane di rections, the flow alig ns with the applied force. A similar conclusion has been 



m-plane dir ections, tne now alig ns with the applied lorcc. A similar conclusion has been 
reached bv lKamrin et al. ( 2010l ). The authors constructed an asymptotic series for 3D 
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flow u over an arbitrary weakly slipping 2D texture with a locally varying slip length 

f3(x,y) : 

OO 

u = J2eW BS - (2-5) 

fe=i 

For ID textures, the boundary conditions to u£ BS can be readily obtained by substitut- 
ing Eq. (|2.5j) into Eq. (|2.3I) and by collecting the terms of the order of e k : 

z = 0: u? BS = f3(x,y), 

k>l: u% BS = p{x,y)d z u™ s . (2 ' 6) 

If we consider periodic textures with j3 (y) being an even function, the slip length can 
be expanded as a cosine Fourier series: 

P= y + ^rcos(27rny), (2.7) 

n=l 
,1/2 

p {y) cos (2irny) dy. (2.8) 

'1/2 

The expansions of the effective slip leng ths for ID textures u p to second order in e of the 
effective slip lengths are then given by ( Kamrin et al. 20101 ) 



^ B /L = e^--e\'£n\p n , (2.9) 



71=1 

6 c J ff A = £Y-e 2 2 7 r^n|/3" . (2.10) 

n=l 

The first-order terms are the isotropic part of the effective slip, Eq. (jl.lj) . The second- 
order terms, which can be neglected for a weakly slipping patterns, are expected to 
introduce the influence of the surface structure, and are responsible for the anisotropy of 
the flow. 

The argument presented above implicitly assumes that the infinite sums over n in the 
higher-order expansion coefficients converge, which implies that the Fourier series, Eq. 
(|2.7p . can be differentiated infinitely often with respect to y. In cases of discontinuous slip, 
where j3(y) exhibits jumps, this is no longer correct and the argument breaks down. From 
a physical point of view, the problem can be associated with singularities of the velocity 
gradient at the boundaries of the slip region. As a specific example, let us consider a 
classical case of alternating weakly slipping (e = bo/L <C 1) stripes with 

B(v)= I 1 aS |y| " 0/2 (2 11) 

P[y > { as 0/2 < \y\ < 1/2 ' 1 > 

so that the boundary conditions, Eq. (|2.3p . can be rewritten as 

z = 0, \y\ < 0/2 : u - ed z u = e, . , 

0/2<|y|<l/2: u = 0. { > 

The Fourier coefficients for this texture then follow from Eq. (|2.8j) 

/3° = 20, (2.13) 
2 sin (7m0) 



n > ; f3 n = 



7m 
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This implies that series in Eqs. (|2.9j) and (|2.10[) . 



71=1 



oo 1 

E 1 



cos (2irnq 



(2.14) 



diverge, since their terms decay as n 1 at n — > oo (large n correspond to small length 



scales). The slow decay of 



with n and the divergence of the scries indicate that 

the expansion (|2.5[) does not resolve properly the solution at small length scales. The 
velocity gradient grows infinitely near the edge of the slip region (see Sec. 13.21 for a 
detailed analysis). As a result, the corresponding term in Eq. (|2.12[) . ed z u, has the same 
order of magnitude as u in the vicinity of the slipping boundary. Therefore, it cannot be 
neglected compared to the leading order, even though e is small. 



3. Weakly slipping stripes 

3.1. Slip-length tensor. 

We now consider the case of stripes more specifically. We first compute the eigenvalues 
of the effective slip-length tensor. Since we assume only weak local slippage, we evaluate 
the effective slip length in the principal directions to second order in e and seek for a 
solution which is finite, i.e., has no singularity. 

A general solution satisfying the Laplace equation (12.21) and decaying at infini ty can 
be presented in terms of a cosine Fourier series as (jAsmolov fc Vinogradovall2012l ) 

a °° 

u = — + 2^ a " CX P (— 27mz) cos (2imy) , (3.1) 

n=i 

where a n are constant coefficients to be found from (|2.12[) . We construct the asymptotic 
series for alternating stripes, 

oo oo 

fe=l fc=l 

imposing that \uk+i/uk\ <C 1 over the entire flow region. The boundary conditions for 
Uk at the wall can be chosen as follows: 

z = : u k - ed z u k = r k (y) , (3.2) 

1 \ as 4>/2 < \y\ < 1/2 ' [ ' 

k>l: r k = l ° ^ aL~ V\<r 1 /o ■ (3-4) 

\ -edzUk-x as (f>/2 < \y\ < 1/2 y J 

The reader may check by the summation of Eqs. (|3.2j) over k that they are fully equivalent 
to Eq. (|2~T2)) . 

The slip velocity is then 

oo 

Usiip = ^a° k /2. 
fc=i 

The boundary condition p.2[) can be rewritten in terms of Fourier coefficients as 



oo 
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The coefficients aj? are now determined using the inverse Fourier transform: 



4 = 2 



n>0 



1/2 
1/2 

2 



rk (y) dy, 

1/2 



(3.5) 



1 + 2-nen 



1/2 



r k (y) cos (2irny) dy. 



From Eqs. 



and (|3.5[) . we have to leading order in e 



n > 



2e4>, 
2e sin (jrnq 



nn (1 + 2-Ksn) 

To find the second-order terms we must evaluate T2 = d z u\, which gives 

sin (Tm(f>) cos (2imy) 



d z u x 



n=l 



a"27rn cos (2irny) = — 4e 



71=1 



1 + 271-en 



The second-order slip velocity is then 

,1/2 



0/2 



8e ■» 
ed z uidy = ^ : 

7T 



7T 



OO 



^ n (1 

n — 1 v 

In (27re) - 



27ren 
1 



sm (nnq) ) 
n (1 + 2iren) 

n—l v 7 

^\ cos (27rn</>) 

n=l 



0(e) 



(3.6) 



(3.7) 



(3.8) 



7' 



In [4 sin 2 (tt^)] | + O (e 3 ) , (3.9) 



where 7 = 0.5772157... is Euler's constant. The series in Q3.8P are very similar to those 
in (|2.14[) . They differ only by the factor (1 + 27ren) in the denominator of the first sum. 
This factor is small at n ~ 1, but it grows linearly with n at large n, thus ensuring 
convergence of the series. Note that the first logarithmic term in (|3.9p does not depend 
on the fraction of the slip regions, <fi. This term is associated with the flow singularities 
near the boundaries between no-slip and slip regions (see Subsection I3.2[) . which are 
responsible for additional viscous dissipation that reduces 6 e ff- 

Finally, for the longitudinal effective slip we obtain the following expansion 



2c 2 

7T 



In 



TVS 



sin (7T<5 



7 



0(e A me) 



from which we can derive the transverse effective slip using (|2.ip . 

2lT£ 



4e 2 



In 



sin (7T<5 



7 



O e d lne 



(3.10) 



(3.11) 



To summarize, we have here directly demonstrated that Eq. (|1.1[) must be applied with 
care. On the one hand, Eqs. (|3.10[) and p. lip unambiguously show that Eq. (jTTTJ) docs 
indeed give the correct first-order term of an expansion for the eigenvalues of the slip- 
length tensor, even in a case of alternating stripes. On the other hand, the higher order 
contributions may be nonanalytical in e, which may create complications. In case of a 
local slip which exhibits step-like jumps at the edge of heterogeneities, the second-order 
terms of the expansions become of the order of e 2 lne (in contrast to e 2 , which would be 
expected for continuously varying local slip). Therefore, they can be comparable to the 
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first-order terms and cannot be ignored even at relatively small e (see Section [5]). These 
terms are not only responsible for anisotropy of the flow, but also (being negative) for 
an additional dissipation. 



3.2. Edge singularity 



We no w describe the flow singularities near slipping heterogeneities in more detail. [Wang 



(|2003h has considered the flow over a surface with rectangular grooves, and found that the 
shear stress is singular near sharp corners, i.e., proportional to r -1 / 3 for longitudinal and 
to r- - 455 ^ for transverse configurations. Here r is the distance from the corner. Following 
this approach, we now consider the flow in the vicinity of the edge of our weakly slipping 
regions, by using polar coordinates (r, 6) with the origin in (y, z) = (<p/2,0). The no-slip 
and slip regions then correspond to 9 = and 9 = tt. The solution of the Laplace equation 
(12. 2j) that satisfies the no-slip boundary condition at 9 — is 

u = cr x sin (A6>) . (3.12) 

The velocity at the edge is finite provided A > 0. The components of velocity gradient 
are 

d z u = cXr x ~ 1 cos [9 (1-X)], d y u = -cAr A ~ 1 sin [9 (1 - A)] . (3.13) 

The velocity decays faster than its gradient as r — > : r x vs. r 1 . Hence, in a small 
region r ~ e, the dimensionlcss shear rate ed z u is of the same order as u, and cannot 
be ignored in the boundary condition (|2.12[) even though e<1. Moreover, at smaller 
distances, r <C e, the term ed z u dominates over it, and the condition in this region 
becomes shear-free: 

r < e, 9 = tt : d z u = 0. (3.14) 

The last condition enables us to find A. To satisfy Eq. (|3.14p one should require, in view 
of p,13p . A = 1/2. Therefore, the velocity over the slip region is 

r<e, 9 = n: u = cr 1/2 , d y u = -cr- 1 ' 2 /2, (3.15) 

where c is a constant. The velocity gradient over the no-slip region follows from (|3.13|) : 

r < e, 9 = : d z u = cr~ 1/2 /2. (3.16) 

In other words, the shear stress has a singularity at the edge. 

We remark that Eqs. (|3.15[) and (|3.16p are valid in a small region, r <C b/L, r <C 1, near 
a jump in the discrete local slip length, from to a finite b. Therefore, our asymptotic 
theory is only valid provided that the fractions of the slip and no-slip regions are not too 
small, 1 — <f> ^> e. Otherwise, the two edges of heterogeneities are close to each 

other, so that the singular regions overlap. Note that a similar, r -1 / 2 , dependence of the 
velocity has been obtained earlier for a no-slip surface decorate d with perfect-slip stripes 
( Sbragaglia fc Prosperetti 2007 ; Asmolov fc Vinogradova 20121 ). 



For the tr ansverse flow one can use the re lation between the velocity fields for the two 



orientations ([Asmolov fc Vinogradovall2012[ ) 



1 / . du d \ z du d 

2{ Ud + Z ^z-)> W = (3J7) 

P = - 9 i±, (3.18) 
dy 



where Ud (y, z) = u [y, z, 2e/3 (y)} is the velocity field for the longitudinal pattern with 
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double local slip length (cf. Eq. (|2.1j) . Hence we conclude that at the wall 

1 dv du d 

2 = : v= 2 Ud > 9i=aT (3 - 19) 

From Eqs. Q3.15P - p,18[) . it also follows that and p all have the same singularity 

r^ 1 / 2 at the edge of weakly slipping region. 

4. Simulation method 

Ou r simulations are performed using Dissipative Pa rticle Dynamics ( Koelman fc Hoogerbruggel 



1993 HEspafiol fc WarrerJ 19951 iGrootfe Warrenll 19971 ). an established method for mesoscale 



fluid simulations. The DPD method is a coarse-grained, momentum-conserving simula- 
tion technique, which naturally includes thermal fluctuations. More spec ifically, we use 



a DPD version without conservati ve interactions (ISoddemann et al\ 120031 ) and combine 



that with a tunable-slip method (jSmiatek et al\ 120081 ) which allows one to implement 



arbitrary hydrodynamic boundary condition. 

The basic DPD equations involve pair interaction between fluid particles. The force 
exerted by particle j on particle i is given by 

F° PD = F% + F*. (4.1) 
The dissipative part is proportional to the relative velocity between two particles, 

F° = -jdpd un{rij)(vij ■ r i:j )r tj (4.2) 

with the friction coefficient ^dpd- The weight function cj£)(r.jj) is a monotonically de- 
creasing function of and vanishes at a given cutoff r c . The random force Ffj has the 
form 



Ffj = y]2k B T'iDPDU D {r ij )^jrij, (4.3) 

where £y = £jj are symmetric, but otherwise uncorrelated random functions with zero 
mean and variance (£ij(t)f;ij(t')) = S(t—t') (here S(t) is Dirac's delta function). The mag- 
nitude of the stochastic contribution is related to the dissipative part by the fluctuation- 
dissipation theorem to ensure correct equilibrium statistics. The pair forces between two 
particles satisfy Newton's third law, = —Fji, and hence the momentum is conserved. 
This leads to correct long-time hydrodynamic behavior (i.e. Navier-Stokcs). 

The wall interaction is introduced in the same spirit. The force on particle i from the 
channel is given by 

F wall = pWCA + p D + F fl ( 4 _ 4 j 

The first one is a repulsive interaction to prevent the fluid particles from penetrating the 
wall. It can be written in terms of the gradient of a Weeks-Chandler- Andersen (WCA) 
potential 

,r»~v.v W , n^p^-p + S (4.5) 

where z is the distance between the fluid particle and the wall. The dissipative contri- 
bution is similar to Eq. (|4.2[) . with the velocity difference v%j replaced by the particle 
velocity relative to the wall, 

Ff = - lL uj L {z)(v t - v waU ). (4.6) 



The parameter characterizes the strength of the wall friction and can be used to tune 
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Fluid density p 3.75a -3 

Friction coefficient for DPD interaction ~/dpd 5.0^/me/cr 

Cutoff for DPD interaction r c 1.0 a 

Coefficient 7l for no-slip boundaries 5.26^/me/a 

Cutoff for wall interaction 2.0 a 

Coefficients jl for slip boundaries {0.5, 0.4, 0.3, 0.2, 

Corresponding slip lengths b {2.6, 3.4, 4.7, 7.1, 

Shear viscosity n s 1.35 ± Q.Qly/me/a 2 

Position of hydrodynamic boundary 1.06 ± 0.12cr 

Table 1. Parameters used in the DPD simulations. 



the value of the slip length. For example, 7l = corresponds to perfectly slippery wall, 
while a positive value of jl leads to a finite slip length. The random term has the form 



Ff = y/2k B T lL Lj L (z)^, (4.7) 

where each component of £ 4 is an independent random variable function with zero mean 
and variance a = — *')• 

The simulations ar e carried out using the open source simulation package ESPResSo 
( Limbach et oi.ll2006t ). We use a quadratically decaying weight function wrj(r) and a 



linearly decaying weight function cjj,(r). Table Q] summarizes the simulation parameters. 
The most important quantity is the slip length b, wh ich can be est i mated from the 
simulation parameters using an analytical formula due to Smiatek et al. ( 2008t ). However, 



the accuracy of the analytic prediction is not satisfactory for weakly slipping surfaces. 
We simulated plane Poiseuille and Couette flows with various 7l to obtain t he slip length 



and the position of the hydrodynamic boundary (see lSmiatek et all (|2008f ) for details). 
Poiseuille flow was implemented by applying a constant body force of the order 10~ 4 e/o 
to all particles. 

In the simulations, the patterned surface has a stripe spacing of L = 100cr. The simula- 
tion box is a rectangular cuboid of size 20a x 100<7 x 102<t. Periodic boundary conditions 
are used in the xy plane. The integration time step is O.Oli/m/ecr. Due to the large sys- 
tem size, the time for reaching a steady state is also quite large (over 10 6 time steps). For 
weakly slipping surfaces, the flow velocity near the wall is small compared to the thermal 
fluctuations; thus long simulation times are also required to obtain enough statistics. In 
this work, velocity profiles have been averaged over 2 x 10 5 time steps. 

Based on the values of the velocities close to the surface, we can estimate the charac- 
teristic Reynolds number of our system. Typical Reynolds numbers in our simulations 
are of order 10, which is much larger than typical values in microfluidic devices. Thus, 
inertia effects may become important in the simulations, and the Stokes equation is not 
strictly valid. We shall sec below that this slightly affects our simulation results. To reach 
more realistic Reynolds numbers, one would have to reduce the shear rate by four orders 
of magnitude. This would reduce the average flow velocity significantly, and the neces- 
sary simulation time to gather data with sufficiently good statistics will then increase 
prohibitively. 
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(a) o.o3o 



0.025 
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Figure 2. The longitudinal effective slip length as functions of (a) the local slip for the texture 
with <j> = 0.3, and (b) the fraction of the slipping phase at b/L — 0.071. Symbols are simulation 
data. Dash-dotted curves show exact numerical results, solid lines correspond to the two-term 
logarithmic expansions (Eq. (|3.10|> ). dash lines to the linear theory (Eq. (|l.ip ). dotted lines to 
Eq. (Onj- 



5. Results and Discussion 

In this section, we compare predictions of our asymptotic theory with results of DPD 
simulations and direct numerical solutions of Eqs. (|2.2[) - (|2.4[) (sec Appendix [S] for details 
of our numerical scheme) . 

Figures HJa) and (b) show the exact numerical results and DPD simulation data for the 
longitudinal component of the slip-length tensor as a function of the dimensionlcss slip 
length b/L and the slipping area <p. The simulation data are in excellent agreement with 
the numerical results, confirming the validity of our DPD scheme. Similar calculations 
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were made for the transverse component of the slip-length tensor. All curves were found 
to be very similar to those presented in Figures [21 therefore, we do not show them here. 
The values for the transverse component are smaller than those for the longitudinal 
component, indicating that the flow is anistropic. The simulation data in the transverse 
case tend to be slightly small er than the predict ion from the numerical solution. This 
has been observed previously IZhou et al (2012) and can be related to the relatively 
large Reynolds numbers in our system (see Sec. [3]). Inertia effect influence the flow past 
transverse stripes, as will be discussed below in the context of Figure 3J For flow past 
longitudinal stripes, the inertia effects are much smaller. Thus the DPD data shown in 
Fig. [2] are not affected by Reynolds number effects. 

The surface-averaged slip, predicted by Eq. (|1.1[) , is also shown in Figure [2] and is well 
above the exact values of the longitudinal effective slip. Also included in Figure[2]are the 
predictions of our theoretical result, Eq. (|3.10|1 . One can see that Eq. (|3.10j) indeed gives 
the correct asymptotic behavior in the limit of very small b/L. It slightly overestimates 
the value of the longitudinal effective slip at larger b/L. We remark and stress that 
nevertheless , our second-order calculation is much more accurate than Eq. (|l.ll) . 

Recently, iBelvaev fc Vinogradoval (|2010al ) suggested approximate expressions for ef- 
fective slip lengths of a surface decorated by partial slip stripes: 
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These formulae have been verified aev fc Vinogradoval l2010al ) using the method 



developed by ICottin-Bizonne et all (|2004l ). The agreement between the theoretical and 
numerical data was found to be very good for all (j> and b/L, but at b/L = O(l), small 
discrepancies were observed, suggesting that Eqs. (|5.1|) and (|5.2[) slightly underestimate 
the effective slip length. To examine this more closely, we also include the prediction of 
Eq. (|5.1|) in Figure [2] We find indeed a small discrepancy between the exact numerical 
data and the predictions of Eq. (|5.1[) . which gives smaller values for the slip length. The 
same trends were observed in a wide range of 0, and the discrepancy slightly increases 
wi th the fraction of slippin g phase . Still, the analytical expressions for the effective slip 
by IBelvaev fe Vinogradova ( 2010a ) appear to be surprisingly accurate, given their sim- 
plicity. We stress, however, that they do not reproduce the asymptotic result, Eq. (II. ip . 
in the limit of very small b/L. They do correctly predict a linear dependence on b in the 
limit of weakly slipping stripes 



but the prefactor, /(</>), differs from 
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This prefactor corresponds to the slope of the curve b oS (b) at b/L = 0. In Figure O 
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Figure 3. The difference between iongitudinai and transverse effective siip lengths as functions 
of (a) the local slip for textures with <j> = 0.3, and (b) <f> for textures with b/L — 0.071. Solid 
curves correspond to calculations made with the two-term logarithmic expansions [Eqs. 1)3.101) 
and (|3.11[) ]. Dotted curves are obtained using Eqs. (|5.1[) and (|5.2[) . Other notations are the same 
as in Figure 



the slope of the dotted line corresponding to Eq. ()5.1|) is smaller than the exact one. 
Nevertheless, the values of &jj ff and b^ s given by (|5 . 1 1) and (|5.2[) correlate well with the 
numerical data for all <p and small but finite b/L. 

It has been predicted that in a thick channel the amplitude of a transverse flow, 
or the flow anisotropy, is controlled by the difference between the eigenvalues of the 
effect ive slip tensor, b\ s — b^r s , which increases with </> and b ( Vinogradova fc Belvaev 



20111 ). According to Eq. (|l.ip . this difference should vanish for weakly slipping surfaces. 



The effect of anisotropy is highlighted in Figure EJa) and (b), which shows the difference 
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between the longitudinal and transverse effective slip lengths computed for fixed <f> = 0.3 
and b/L = 0.071, respectively. The exact numerical values are positive, except in the 
case of extremely small local slip, clearly showing that the flow is anisotropic. This 
is confirmed by the simulation results. The error bars are relatively large. For weakly 
slipping surfaces, the difference &j| ff — b^ s is small compared to the slip lengths themselves, 
(of the order of e 2 lne), and this is the reason for the large error of the simulation data. 
The simulation data agree with the numerical results within the error. Nevertheless, the 
data suggest that they lie systematically above the numerical results especially for larger 
slipping phase fraction (p. This is a consequence of the relatively large Reynold number. 
As discussed above, inertia effects primarily affect the flow and effective slip length in 
the transverse configuration. Test runs with larger shear rates were performed, and the 
deviations increased, indicating that they presumably vanish in the Stokes limit. 

The results from the two-term logarithmic expansions, Eqs. (|3.10[) and (|3.11[) , are also 
shown in Figure [3] In the limit of small b/L, the two-term expansion predicts correctly the 
positive difference and enhanced anisotropy as the slip length increases. At larger b/L, 
deviations from the numerical results become larger due to the increasing contribution 
from higher order terms. At b/L = 0.071, the two-term prediction for the slip length 
difference is only in moderately good agreement with the numerical data. For very low or 
very high coverage (4> — > or (f> — > 1), the agreement is not good at all, the theory even 
predicts the wrong sign [Figure [3jb)]. This is consistent with our discussion in Section 
13.21 where we have argued that the approximation must break down when the singular 
regions associated with adjacent edges overlap. Also included in Figure [3] is the result 
from the approximate expressions Eqs. (|5.1[) and (|5.2p . which again shows surprisingly 
good agreement with the numerical data over the whole range of (j>. 

Figure SJa) presents the longitudinal velocity at the wall for various b/L. In the sim- 
ulations, the slip velocity has been obtained from an extrapolation procedure. Due to 
the small magnitude of the slip velocity in comparison to the thermal fluctuation (order 
of 1 for ksT = 1), the data scatter very much. Much longer averaging times would be 
necessary to improve the statistics. The agreement of the exact numerical results and 
DPD simulation data is again very good. The velocity distribution is not smooth at the 
edge, instead it rises according to a power law on the slipping area with exponent close 
to A = 0.5 as predicted in Section T3.2I In Figure HJb), we verify the relation between 
the transverse and longitudinal velocities. Eq. (|3.19[) suggests that the local slip veloc- 
ity above transverse stripes of slip length b should be identical to half of that above 
longitudinal stripes of slip length 26. This is confirmed by the numerical results. The 
simulation data, however, show deviations near the edge. This illustrates the origin of 
the finite Reynolds number effects discussed above. In simulations, the fluid is modelled 
as DPD particles with finite mass, and abrupt changes of the transverse velocity are 
suppressed because of inertia. Therefore the transverse velocity is smoothed out near the 
edge, showing a smaller value compared to the numerical data. 



6. Concluding Remarks 

In conclusion, we have investigated shear flow past weakly slipping super-hydrophobic 
stripes, focussing in particular on edge effects associated with steplike discontinuities in 
the local slip length. The essential conclusion from our analysis is that such step ef- 
fects reduce the effective slip below the surface-averaged value and induce anisotropy. In 
practice, this means that the flow does not align with the applied shear str ess. Thus, it 



should be possible to generate transver se hydrodynamic phenomena (like in IStone et 



( 2004 ): Vinogradova fc Belvaev ( 201ll )) even with such weakly slipping anisotropic tex- 
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Figure 4. (a) The longitudinal velocity along the wall for a texture with 4> = 0.7 and differ- 
ent slip lengths: b/L — 0.145 (circles and solid line), b/L = 0.071 (squares and dashed line), 
and b/L — 0.047 (triangles and dash-dotted line). Symbols are simulation data and lines are 
numerical results, (b) Comparison of the transverse velocity (v[y,z = 0, /3(y)]) and longitudinal 
velocity for nearly double local slip length (^Ud[y, z = 0,2f3(y)], see Eq. (|3 19|) \ The simulation 
and numerical results for longitudinal stripes of b/L = 0.147 are shown in circles and solid line, 
respectively. Squares and dashed line correspond to transverse stripes of b/L = 0.071. 



tures . This may also have relevance for transverse e l ectrokinetic phenomen a (Bahg a et al 



2010t I Vinogradova fc Belvaevll201lt ISauiresI [200l iBelvaev fc Vinogradova! l201ll ). Asa 

side remark, our analytical result opens the possibility of solving analytically many fun- 
damental problems involving, weakly slipping heterogeneous surfaces, inc luding hydrody- 
namic interactions ( Belvaev fc Vinogradova 20106 ; Asmolov et adlioilh . 

Finally, we note, that even though our discussion has been limited to weakly slipping 
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heterogeneities, our model is much more general. Every result in this work could be 
used for describing weakly rough or porous surfaces since at large distances from the 
wall, the boundary condition at the rough interface o r fluid-porous interface may be 
approximated by a, slip model (Beavers fc Joseph! 19671: [Taylor 1971 ; Richardson 1971 



Kunert fc Hartindl2007t iKunert et al.ll20ld: iLecoa et a/.ll2004D . 



Appendix A. Numerical method 

The Navier slip boundary condition (|2.3[) can be written in terms of the Fourier coef- 
ficients a n , accounting for (|3.ip . as 

o 00 

— + I 1 + 2nn£ P (»)] «" cos ( 2 ™y ) = £ P (y) ■ ( A 1 ) 

n=l 

To find a n numerically we truncate the sum at some cut-off number N (usually TV = 512) 
and apply a Fourier cosine transform to (|A lj) . The Fourier coefficients A™ and (3 m of the 
vectors [1 + 27rne/3 (yi)] cos {2imyi) and f3 (yi) , y\ = 1/2N, n, m, I = 0, N, are evaluated 
using the IMSL routine FCOST. Then (jXTj) is reduced to a linear system A™ a n = ]i m , 
which is solved using the routine LSARG. In view of (|3.1j) . the application of FCOST to 
a vector a n exp {—2imz) produces an inverse Fourier transform, i.e., we evaluate u (yi, z) 
for given z. 
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